source("~/jules_ppe_gmd/docs/JULES-ES-1p0-common-packages.R")
## Loading required package: spam
## Loading required package: dotCall64
## Loading required package: grid
## Spam version 2.7-0 (2021-06-25) is loaded.
## Type 'help( Spam)' or 'demo( spam)' for a short introduction 
## and overview of this package.
## Help for individual functions is also obtained by adding the
## suffix '.spam' to the function name, e.g. 'help( chol.spam)'.
## 
## Attaching package: 'spam'
## The following objects are masked from 'package:base':
## 
##     backsolve, forwardsolve
## Loading required package: viridis
## Loading required package: viridisLite
## See https://github.com/NCAR/Fields for
##  an extensive vignette, other supplements and source code
## Loading required package: lhs
## Loading required package: parallel
## Loading required package: viztools
source("~/jules_ppe_gmd/docs/JULES-ES-1p0-common-functions.R")


## ----------------------------------------------------------------------
## Data locations and constants
## ----------------------------------------------------------------------
#ensloc <- '/project/carbon_ppe/JULES-ES-1p0_PPE/'
ensloc_wave00 <- '/data/users/hadaw/JULES_ES_PPE/u-au932/'
ensloc_wave01 <- '/data/users/hadaw/JULES_ES_PPE/u-ck006/'

 
# Some pallete options
yg = brewer.pal(9, "YlGn")
ryb = brewer.pal(11, "RdYlBu")
byr = rev(ryb)
rb = brewer.pal(11, "RdBu")
br = rev(rb)
blues = brewer.pal(9, 'Blues')
cbPal <- c("#000000", "#E69F00", "#56B4E9", "#009E73", "#F0E442", "#0072B2", "#D55E00", "#CC79A7")

wave00col <- 'skyblue2'
wave01col <- 'tomato2'

zissou5 <- wes_palette('Zissou1', 5, type = c('discrete', 'continuous'))
zblue <- makeTransparent(as.character(zissou5)[1], 150)
zred <- makeTransparent(as.character(zissou5)[5], 150)

ysec = 60*60*24*365
years <- 1850:2013
global_mean_modern_value_wave00_file <- "~/jules_ppe_gmd/docs/data/global_mean_modern_value_wave00_2022-13-09.rdata"
load(global_mean_modern_value_wave00_file)

modern_value_stan_file <- "~/jules_ppe_gmd/docs/data/modern_value_stan_2022-09-13.rdata"
load(modern_value_stan_file)
  
Y <- datmat_wave00

Y_nlevel0_ix <- which(is.na(datmat_wave00[,'year']))
# Load up the data
lhs_i = read.table('~/jules_ppe_gmd/docs/data/lhs_u-ao732.txt', header = TRUE)
lhs_ii = read.table('~/jules_ppe_gmd/docs/data/lhs_u-ao732a.txt', header = TRUE)

toplevel_ix = 1:499

# The raw input data is a latin hypercube
lhs = rbind(lhs_i, lhs_ii)[toplevel_ix, ]
lhs_level0 <- lhs[-Y_nlevel0_ix,]

X = normalize(lhs)
colnames(X) = colnames(lhs)

X_level0 <- X[-Y_nlevel0_ix,]
X_nlevel0 <- X[Y_nlevel0_ix,]

d = ncol(X)
# lower and higher bound on the normalised matrix for visualisation
rx = rbind(rep(0,32), rep(1,32))
low_npp_ix <- which(Y[,'npp_nlim_lnd_sum'] < 1e5)
# code from https://stackoverflow.com/questions/28182872/how-to-use-different-sets-of-data-in-lower-and-upper-panel-of-pairs-function-in


#X <- matrix(runif(300), ncol=3)
#Y <- matrix(c(sort(runif(100, 0, 10)), 
#              sort(runif(100, 0, 10)), 
#              sort(runif(100, 0, 10))), ncol=3)

#pdf(file = 'figs/fig02.pdf', width = 12, height = 10)
#pdf(file = 'figs/run-failure-pairs.pdf', width = 12, height = 10)
x1 <- X[low_npp_ix, ]
x2 <- X_nlevel0

XY <- rbind(x1, x2)


pairs(XY,
      lower.panel=function(x, y, ...) {
        Xx <- x[seq_len(nrow(x1))] # corresponds to X subset
        Xy <- y[seq_len(nrow(x1))] # corresponds to X subset
        #usr <- par("usr"); on.exit(par(usr))
        #par(usr = c(range(x1[, -ncol(x1)]), range(x1[, -1]))) # set up limits
        points(Xx, Xy, col = zblue, pch = 19, cex = 0.8)
       # if(par('mfg')[2] == 1) axis(2) # if left plot, add left axis
        #if(par('mfg')[1] == ncol(x1)) axis(1) # if bottom plot add bottom axis
      }, 
      upper.panel=function(x, y, ...) {
        Yx <- x[(nrow(x1) + 1):length(x)] # Y subset
        Yy <- y[(nrow(x1) + 1):length(y)] # Y subset
        
        #cntr <- outer(Yx, Yx, FUN='*') # arbitrary function for contour
       # usr <- par("usr"); on.exit(par(usr))
        #par(usr = c(range(x2[, -1]), range(x2[, -ncol(x2)]))) # set up limits
        points(Yx, Yy, col = zred, pch = 19, cex = 0.8)
        #contour(Yx, Yy, cntr, add=TRUE)
        #if(par('mfg')[2] == ncol(x2)) axis(4) # if right plot, add right axis
        #if(par('mfg')[1] == 1) axis(3) # if top plot, add top axis
      }, 
      #tick=FALSE, # suppress the default tick marks
      #line=1,
      gap = 0,
      xlim = c(0,1), ylim = c(0,1),
      labels = 1:d,
      oma = c(2, 18, 2, 2)) # move the default tick labels off the plot 

reset()

legend('left', legend = paste(1:d, colnames(lhs)), cex = 1.1, bty = 'n')
legend('topleft', pch = 19, col = c( zred, zblue), legend = c('failed', 'zero carbon cycle'), bty = 'n', inset = 0.02, cex = 1.1 )

#dev.off()

Generate the data to save

Y_char <- rep('RAN', nrow(X))

Y_char[Y_nlevel0_ix] <- 'CRASHED'
Y_char[low_npp_ix] <- 'LOWCARBON'

Y_factor <- as.factor(Y_char)

wave0_summary_df <- cbind(X, Y, Y_char, Y_factor)
save(lhs, X, Y, Y_char, Y_factor, wave0_summary_df,  file = 'data/wave0_summary.RData')

Data checks

Have we got the factors right?

cbind(Y[, 'npp_nlim_lnd_sum'], Y_char)
##                             Y_char     
##   [1,] "2464953.47849605"   "RAN"      
##   [2,] "1248000.44246179"   "RAN"      
##   [3,] "1893850.67745132"   "RAN"      
##   [4,] "1779634.61879458"   "RAN"      
##   [5,] "1333052.53686655"   "RAN"      
##   [6,] "1829087.04066031"   "RAN"      
##   [7,] "125.100277029329"   "LOWCARBON"
##   [8,] "1354509.40828833"   "RAN"      
##   [9,] "1436905.38963356"   "RAN"      
##  [10,] "1105586.84200311"   "RAN"      
##  [11,] "2281524.08799275"   "RAN"      
##  [12,] "1569806.15413172"   "RAN"      
##  [13,] "1452198.55278638"   "RAN"      
##  [14,] "1199097.02939196"   "RAN"      
##  [15,] "2617350.39816265"   "RAN"      
##  [16,] "1249447.10148417"   "RAN"      
##  [17,] "939242.694289745"   "RAN"      
##  [18,] "5.37100446251452"   "LOWCARBON"
##  [19,] "1077579.05514713"   "RAN"      
##  [20,] "1523180.02583773"   "RAN"      
##  [21,] "1569417.63070336"   "RAN"      
##  [22,] "1873231.77735771"   "RAN"      
##  [23,] "1912583.54670419"   "RAN"      
##  [24,] "1902959.37750349"   "RAN"      
##  [25,] "1445375.9153619"    "RAN"      
##  [26,] "951627.465373486"   "RAN"      
##  [27,] "1662356.03945454"   "RAN"      
##  [28,] "1936958.91626488"   "RAN"      
##  [29,] "1433603.10792285"   "RAN"      
##  [30,] "976651.168260095"   "RAN"      
##  [31,] "2082565.02331613"   "RAN"      
##  [32,] "1531283.77037663"   "RAN"      
##  [33,] "1645926.53544732"   "RAN"      
##  [34,] "1320236.28995648"   "RAN"      
##  [35,] "2891640.62239662"   "RAN"      
##  [36,] "1635834.14721978"   "RAN"      
##  [37,] "2089155.84530289"   "RAN"      
##  [38,] "1377545.60029833"   "RAN"      
##  [39,] "1487813.03878515"   "RAN"      
##  [40,] "651768.709997753"   "RAN"      
##  [41,] "2027624.59401214"   "RAN"      
##  [42,] "2078839.96826463"   "RAN"      
##  [43,] "518555.798220461"   "RAN"      
##  [44,] "1061709.94504286"   "RAN"      
##  [45,] "1775343.21596816"   "RAN"      
##  [46,] "2311775.26816144"   "RAN"      
##  [47,] "1441874.46802137"   "RAN"      
##  [48,] "2195581.99198865"   "RAN"      
##  [49,] "1325496.4911125"    "RAN"      
##  [50,] "2004112.56580796"   "RAN"      
##  [51,] "2368885.55711663"   "RAN"      
##  [52,] "1361.54882660815"   "LOWCARBON"
##  [53,] "908731.084399189"   "RAN"      
##  [54,] "1442169.69546151"   "RAN"      
##  [55,] "139.629393179477"   "LOWCARBON"
##  [56,] "1333927.55011543"   "RAN"      
##  [57,] "1666849.46116248"   "RAN"      
##  [58,] "965867.767068223"   "RAN"      
##  [59,] "1243866.85654177"   "RAN"      
##  [60,] "1943266.73141755"   "RAN"      
##  [61,] "1902138.05825955"   "RAN"      
##  [62,] "1767409.94121306"   "RAN"      
##  [63,] "1317749.66685722"   "RAN"      
##  [64,] "1199841.24949906"   "RAN"      
##  [65,] "1837154.16360185"   "RAN"      
##  [66,] "1163145.1259366"    "RAN"      
##  [67,] "1214317.39561894"   "RAN"      
##  [68,] "1313614.75168373"   "RAN"      
##  [69,] "1724789.67747131"   "RAN"      
##  [70,] "1576737.37471207"   "RAN"      
##  [71,] "781856.046829388"   "RAN"      
##  [72,] "1101375.98446798"   "RAN"      
##  [73,] "1817092.60446859"   "RAN"      
##  [74,] "2105.24224905522"   "LOWCARBON"
##  [75,] "1826534.74527454"   "RAN"      
##  [76,] "1155480.88628367"   "RAN"      
##  [77,] "1471672.36197739"   "RAN"      
##  [78,] "633901.745067233"   "RAN"      
##  [79,] "911610.216606713"   "RAN"      
##  [80,] "2290213.70950657"   "RAN"      
##  [81,] "2646262.40850179"   "RAN"      
##  [82,] "1352243.52082662"   "RAN"      
##  [83,] "1265953.96448741"   "RAN"      
##  [84,] "2100616.21162727"   "RAN"      
##  [85,] "1776678.07569972"   "RAN"      
##  [86,] "836.397054014272"   "LOWCARBON"
##  [87,] NA                   "CRASHED"  
##  [88,] "438550.74712463"    "RAN"      
##  [89,] "1279.62367587908"   "LOWCARBON"
##  [90,] "1524159.38243703"   "RAN"      
##  [91,] "2188559.27733915"   "RAN"      
##  [92,] "1050366.8900468"    "RAN"      
##  [93,] "1429623.48469536"   "RAN"      
##  [94,] "27456.120078662"    "LOWCARBON"
##  [95,] "1398.90443438924"   "LOWCARBON"
##  [96,] "1037573.57180142"   "RAN"      
##  [97,] "80867.6892480567"   "LOWCARBON"
##  [98,] "1744445.9170165"    "RAN"      
##  [99,] "1233881.05138578"   "RAN"      
## [100,] "1443212.18749931"   "RAN"      
## [101,] "1898806.28032103"   "RAN"      
## [102,] "22.6523656647353"   "LOWCARBON"
## [103,] "1618940.96335165"   "RAN"      
## [104,] "1773612.19063594"   "RAN"      
## [105,] "1505628.99582318"   "RAN"      
## [106,] "1511165.57724497"   "RAN"      
## [107,] "2387524.55902676"   "RAN"      
## [108,] "1808140.75254821"   "RAN"      
## [109,] "3034612.74359176"   "RAN"      
## [110,] "1305467.18461631"   "RAN"      
## [111,] "2215.47450906594"   "LOWCARBON"
## [112,] "339510.74735257"    "RAN"      
## [113,] "1073457.42123026"   "RAN"      
## [114,] "1004856.00646112"   "RAN"      
## [115,] "1580083.43485615"   "RAN"      
## [116,] "921086.751260579"   "RAN"      
## [117,] "897774.501909049"   "RAN"      
## [118,] "2678468.32019488"   "RAN"      
## [119,] "975974.054951735"   "RAN"      
## [120,] "1304973.40986444"   "RAN"      
## [121,] "2795729.48708828"   "RAN"      
## [122,] "972523.559775463"   "RAN"      
## [123,] "1049919.96813186"   "RAN"      
## [124,] "2997.38035283134"   "LOWCARBON"
## [125,] "561.565402536467"   "LOWCARBON"
## [126,] "2988253.78284404"   "RAN"      
## [127,] "2081895.14453078"   "RAN"      
## [128,] "3070797.46549112"   "RAN"      
## [129,] "1403506.15570796"   "RAN"      
## [130,] "827048.517335001"   "RAN"      
## [131,] "395143.022994698"   "RAN"      
## [132,] "2163311.88964093"   "RAN"      
## [133,] "2178050.99534827"   "RAN"      
## [134,] "566076.622930143"   "RAN"      
## [135,] "1392961.37882784"   "RAN"      
## [136,] "18392.8206574553"   "LOWCARBON"
## [137,] "1947440.97630319"   "RAN"      
## [138,] "2636818.64751266"   "RAN"      
## [139,] "28499.059770153"    "LOWCARBON"
## [140,] "1164715.18341187"   "RAN"      
## [141,] "0"                  "LOWCARBON"
## [142,] "835144.438570216"   "RAN"      
## [143,] "2128420.55576015"   "RAN"      
## [144,] "1521308.9730314"    "RAN"      
## [145,] NA                   "CRASHED"  
## [146,] "2748024.4267935"    "RAN"      
## [147,] "973221.738622719"   "RAN"      
## [148,] "860049.23757151"    "RAN"      
## [149,] "1324456.75170062"   "RAN"      
## [150,] "1201010.03107696"   "RAN"      
## [151,] "1142603.5225206"    "RAN"      
## [152,] NA                   "CRASHED"  
## [153,] "2693611.00987939"   "RAN"      
## [154,] "1069541.52591967"   "RAN"      
## [155,] "2393741.91555386"   "RAN"      
## [156,] "1382026.59899902"   "RAN"      
## [157,] "1862826.09944603"   "RAN"      
## [158,] "819708.358546789"   "RAN"      
## [159,] "16.1599645777015"   "LOWCARBON"
## [160,] "10.6491420068659"   "LOWCARBON"
## [161,] "1237522.03995008"   "RAN"      
## [162,] "851733.274253283"   "RAN"      
## [163,] "1499290.2664036"    "RAN"      
## [164,] "969437.761077347"   "RAN"      
## [165,] "654971.149328721"   "RAN"      
## [166,] "1182209.44900622"   "RAN"      
## [167,] "2118636.567432"     "RAN"      
## [168,] "1990818.87407025"   "RAN"      
## [169,] "2257432.14856896"   "RAN"      
## [170,] "1550902.99283446"   "RAN"      
## [171,] NA                   "CRASHED"  
## [172,] "8616.48535791419"   "LOWCARBON"
## [173,] "2063335.56382181"   "RAN"      
## [174,] "2217.33780499236"   "LOWCARBON"
## [175,] "1455273.15033427"   "RAN"      
## [176,] "1602432.8716883"    "RAN"      
## [177,] "1442763.37340397"   "RAN"      
## [178,] "935002.418487836"   "RAN"      
## [179,] NA                   "CRASHED"  
## [180,] "2398227.96374005"   "RAN"      
## [181,] "2257472.96108099"   "RAN"      
## [182,] "1463089.42164426"   "RAN"      
## [183,] "968805.302928887"   "RAN"      
## [184,] "724344.993211166"   "RAN"      
## [185,] "202447.955948166"   "RAN"      
## [186,] "1436879.80366549"   "RAN"      
## [187,] "1050598.86322165"   "RAN"      
## [188,] "858303.300545437"   "RAN"      
## [189,] "837.8470560364"     "LOWCARBON"
## [190,] "77.7376830697535"   "LOWCARBON"
## [191,] "2666235.23760768"   "RAN"      
## [192,] "1442330.64357137"   "RAN"      
## [193,] "300301.613025114"   "RAN"      
## [194,] "1867448.54134229"   "RAN"      
## [195,] "2318123.19343285"   "RAN"      
## [196,] "1192558.07059905"   "RAN"      
## [197,] "296580.206341884"   "RAN"      
## [198,] "198948.53460339"    "RAN"      
## [199,] "1475476.98479445"   "RAN"      
## [200,] "1151974.92734534"   "RAN"      
## [201,] "1045780.41572662"   "RAN"      
## [202,] "597.349982204987"   "LOWCARBON"
## [203,] "2088874.69446935"   "RAN"      
## [204,] "1001228.25695345"   "RAN"      
## [205,] "1163274.13917084"   "RAN"      
## [206,] NA                   "CRASHED"  
## [207,] "16.9065379434552"   "LOWCARBON"
## [208,] "1719257.89125913"   "RAN"      
## [209,] "1223443.31154094"   "RAN"      
## [210,] "436285.944826903"   "RAN"      
## [211,] "2517584.52687686"   "RAN"      
## [212,] "1147738.92352667"   "RAN"      
## [213,] "16630.6275934419"   "LOWCARBON"
## [214,] "1136523.58134658"   "RAN"      
## [215,] "2945929.85082532"   "RAN"      
## [216,] "1027358.29145306"   "RAN"      
## [217,] "1090954.7293875"    "RAN"      
## [218,] NA                   "CRASHED"  
## [219,] "474.410621824389"   "LOWCARBON"
## [220,] "1182142.21931005"   "RAN"      
## [221,] "1244732.72703289"   "RAN"      
## [222,] "817223.408383973"   "RAN"      
## [223,] "789494.766979139"   "RAN"      
## [224,] "1286529.489797"     "RAN"      
## [225,] "771006.894030684"   "RAN"      
## [226,] "778690.42112912"    "RAN"      
## [227,] "177.882608699528"   "LOWCARBON"
## [228,] "2084474.06576161"   "RAN"      
## [229,] "1225.95556037021"   "LOWCARBON"
## [230,] NA                   "CRASHED"  
## [231,] "1756439.8256892"    "RAN"      
## [232,] "14878.4358912422"   "LOWCARBON"
## [233,] "309973.949628558"   "RAN"      
## [234,] "8764.0532200681"    "LOWCARBON"
## [235,] "267339.183244262"   "RAN"      
## [236,] "5418.55514931159"   "LOWCARBON"
## [237,] "1441975.92917763"   "RAN"      
## [238,] "709437.349215471"   "RAN"      
## [239,] NA                   "CRASHED"  
## [240,] "259956.693826904"   "RAN"      
## [241,] "88222.8675543907"   "LOWCARBON"
## [242,] "4766.6621059493"    "LOWCARBON"
## [243,] "143.9677977123"     "LOWCARBON"
## [244,] "1137573.83629277"   "RAN"      
## [245,] "655601.191706854"   "RAN"      
## [246,] "594181.861682584"   "RAN"      
## [247,] "488056.40357227"    "RAN"      
## [248,] "1128079.61589957"   "RAN"      
## [249,] "2897145.33655448"   "RAN"      
## [250,] "2726652.51399176"   "RAN"      
## [251,] "21571.315698287"    "LOWCARBON"
## [252,] NA                   "CRASHED"  
## [253,] "1961574.09889377"   "RAN"      
## [254,] "47946.6186543254"   "LOWCARBON"
## [255,] "27076.4135515304"   "LOWCARBON"
## [256,] "2251606.48396996"   "RAN"      
## [257,] "3273070.39797162"   "RAN"      
## [258,] "645886.579584998"   "RAN"      
## [259,] "586088.226578541"   "RAN"      
## [260,] "1156781.11188133"   "RAN"      
## [261,] NA                   "CRASHED"  
## [262,] "1031063.75077349"   "RAN"      
## [263,] NA                   "CRASHED"  
## [264,] "1917263.39494636"   "RAN"      
## [265,] "1615558.13924909"   "RAN"      
## [266,] "2533649.88926883"   "RAN"      
## [267,] "1531972.92447036"   "RAN"      
## [268,] "1509.83741632229"   "LOWCARBON"
## [269,] "2231990.01459277"   "RAN"      
## [270,] "39017.5736788715"   "LOWCARBON"
## [271,] "0.0789806248899001" "LOWCARBON"
## [272,] "2324156.32963896"   "RAN"      
## [273,] "21.1635773171758"   "LOWCARBON"
## [274,] "588731.278489821"   "RAN"      
## [275,] "8434.87646631176"   "LOWCARBON"
## [276,] "1052314.17562924"   "RAN"      
## [277,] "233938.318100738"   "RAN"      
## [278,] "1353008.97324927"   "RAN"      
## [279,] "7.79517493656395"   "LOWCARBON"
## [280,] "2224267.59887657"   "RAN"      
## [281,] "2480640.44093978"   "RAN"      
## [282,] "8334.02080434849"   "LOWCARBON"
## [283,] "0"                  "LOWCARBON"
## [284,] "2112123.26551258"   "RAN"      
## [285,] "3575717.57030972"   "RAN"      
## [286,] "1137829.11665267"   "RAN"      
## [287,] "768796.012526191"   "RAN"      
## [288,] "182756.792189772"   "RAN"      
## [289,] "2807794.29937303"   "RAN"      
## [290,] "26186.6324336018"   "LOWCARBON"
## [291,] "5.16383244612062"   "LOWCARBON"
## [292,] "5272823.4048439"    "RAN"      
## [293,] "3378889.69905921"   "RAN"      
## [294,] "967239.187151969"   "RAN"      
## [295,] "1513055.5460534"    "RAN"      
## [296,] "853512.111943073"   "RAN"      
## [297,] "1195395.23093059"   "RAN"      
## [298,] "1191767.22530273"   "RAN"      
## [299,] "2981.94763769056"   "LOWCARBON"
## [300,] "1564848.16964212"   "RAN"      
## [301,] "2177454.63091984"   "RAN"      
## [302,] "1585858.82741674"   "RAN"      
## [303,] "1402252.58495003"   "RAN"      
## [304,] "1875165.36531898"   "RAN"      
## [305,] "1428993.72501249"   "RAN"      
## [306,] "2041031.91696762"   "RAN"      
## [307,] "1812481.75362187"   "RAN"      
## [308,] "1499581.42629162"   "RAN"      
## [309,] "1398649.12948884"   "RAN"      
## [310,] "1622463.15225339"   "RAN"      
## [311,] "2021446.30744064"   "RAN"      
## [312,] "1528131.98863308"   "RAN"      
## [313,] "1346977.76387442"   "RAN"      
## [314,] "1803078.20336496"   "RAN"      
## [315,] NA                   "CRASHED"  
## [316,] "3600.55093489759"   "LOWCARBON"
## [317,] "2735351.14292133"   "RAN"      
## [318,] "1875802.01944773"   "RAN"      
## [319,] "1011538.95542792"   "RAN"      
## [320,] "1571171.02786292"   "RAN"      
## [321,] "1441178.46552292"   "RAN"      
## [322,] "1449150.98896031"   "RAN"      
## [323,] "634697.85860513"    "RAN"      
## [324,] "1559508.88789955"   "RAN"      
## [325,] "2079199.4549005"    "RAN"      
## [326,] "31732.7657074611"   "LOWCARBON"
## [327,] "44021.6825703995"   "LOWCARBON"
## [328,] "1107846.07986755"   "RAN"      
## [329,] "1622331.38157197"   "RAN"      
## [330,] "2030849.88858994"   "RAN"      
## [331,] "468054.403308293"   "RAN"      
## [332,] "1203175.42133573"   "RAN"      
## [333,] "53313.7875880245"   "LOWCARBON"
## [334,] "1091233.61100635"   "RAN"      
## [335,] "1731524.89117561"   "RAN"      
## [336,] NA                   "CRASHED"  
## [337,] "2281894.39736456"   "RAN"      
## [338,] "1856521.22642198"   "RAN"      
## [339,] "21.6520858369119"   "LOWCARBON"
## [340,] NA                   "CRASHED"  
## [341,] "2285885.08156568"   "RAN"      
## [342,] "694644.015393941"   "RAN"      
## [343,] "2571.46567734245"   "LOWCARBON"
## [344,] NA                   "CRASHED"  
## [345,] "1864233.90728213"   "RAN"      
## [346,] "1759935.59114204"   "RAN"      
## [347,] "1119324.48885291"   "RAN"      
## [348,] "2252758.98892193"   "RAN"      
## [349,] "1625761.18856225"   "RAN"      
## [350,] "781589.125235808"   "RAN"      
## [351,] "1244575.03675289"   "RAN"      
## [352,] "2355091.83872167"   "RAN"      
## [353,] "675511.888141798"   "RAN"      
## [354,] "1098030.76399517"   "RAN"      
## [355,] "1273025.56057886"   "RAN"      
## [356,] "836188.644024318"   "RAN"      
## [357,] "2597925.33982223"   "RAN"      
## [358,] NA                   "CRASHED"  
## [359,] "679816.768584808"   "RAN"      
## [360,] "1211960.97253689"   "RAN"      
## [361,] "143589.482061787"   "RAN"      
## [362,] "1087122.24802111"   "RAN"      
## [363,] "2382389.99621423"   "RAN"      
## [364,] "1809180.07632715"   "RAN"      
## [365,] "2249280.13648273"   "RAN"      
## [366,] "1882529.41130647"   "RAN"      
## [367,] "5328.34619164727"   "LOWCARBON"
## [368,] "1208986.32211504"   "RAN"      
## [369,] "1107490.30053689"   "RAN"      
## [370,] "2216991.27203092"   "RAN"      
## [371,] "893859.502962013"   "RAN"      
## [372,] "76.2542326522061"   "LOWCARBON"
## [373,] "1312674.28606182"   "RAN"      
## [374,] "415447.916753839"   "RAN"      
## [375,] "850881.513824418"   "RAN"      
## [376,] "885339.13594881"    "RAN"      
## [377,] "2113588.20086514"   "RAN"      
## [378,] "1425799.60275817"   "RAN"      
## [379,] "1128105.10960357"   "RAN"      
## [380,] "2371649.3296751"    "RAN"      
## [381,] "255886.390305109"   "RAN"      
## [382,] "2466945.24454644"   "RAN"      
## [383,] "2033112.04297303"   "RAN"      
## [384,] "1190784.9478732"    "RAN"      
## [385,] "206780.13645465"    "RAN"      
## [386,] "431400.751720456"   "RAN"      
## [387,] "1808897.43954395"   "RAN"      
## [388,] "221742.90645669"    "RAN"      
## [389,] "1499767.91545002"   "RAN"      
## [390,] "825569.031312678"   "RAN"      
## [391,] "2219524.63248659"   "RAN"      
## [392,] "1751328.0315877"    "RAN"      
## [393,] "141464.11415107"    "RAN"      
## [394,] "1905921.73027521"   "RAN"      
## [395,] NA                   "CRASHED"  
## [396,] "86129.1602837143"   "LOWCARBON"
## [397,] "1738959.2987119"    "RAN"      
## [398,] "2400853.04029598"   "RAN"      
## [399,] "595.41911312446"    "LOWCARBON"
## [400,] NA                   "CRASHED"  
## [401,] "2431040.53251091"   "RAN"      
## [402,] "1616335.42542376"   "RAN"      
## [403,] "1412281.70818034"   "RAN"      
## [404,] "1930697.97345133"   "RAN"      
## [405,] "2288.77902682165"   "LOWCARBON"
## [406,] NA                   "CRASHED"  
## [407,] NA                   "CRASHED"  
## [408,] "5679.83400627306"   "LOWCARBON"
## [409,] "2255483.7662744"    "RAN"      
## [410,] "1697200.48172835"   "RAN"      
## [411,] NA                   "CRASHED"  
## [412,] "86772.1408945445"   "LOWCARBON"
## [413,] "58.8614698072582"   "LOWCARBON"
## [414,] "8079.09621326771"   "LOWCARBON"
## [415,] "1222547.82839755"   "RAN"      
## [416,] "1643981.26262162"   "RAN"      
## [417,] "1090488.74884965"   "RAN"      
## [418,] "1844070.58227165"   "RAN"      
## [419,] "1435918.61989897"   "RAN"      
## [420,] "2331558.94159436"   "RAN"      
## [421,] "352790.572217638"   "RAN"      
## [422,] "2612299.38776824"   "RAN"      
## [423,] "2289592.04066901"   "RAN"      
## [424,] "966660.438227921"   "RAN"      
## [425,] "608747.910345104"   "RAN"      
## [426,] "3814774.23854009"   "RAN"      
## [427,] "25651.4463067504"   "LOWCARBON"
## [428,] "859941.950259722"   "RAN"      
## [429,] "1175141.7966137"    "RAN"      
## [430,] "1388753.39016613"   "RAN"      
## [431,] "1307268.4862548"    "RAN"      
## [432,] "1230418.86711069"   "RAN"      
## [433,] "27562.79145666"     "LOWCARBON"
## [434,] "664304.835670389"   "RAN"      
## [435,] "2309620.06944286"   "RAN"      
## [436,] "1442045.89806"      "RAN"      
## [437,] "1080196.27863663"   "RAN"      
## [438,] "12971.8130503748"   "LOWCARBON"
## [439,] "2112311.00279278"   "RAN"      
## [440,] "1865902.63590031"   "RAN"      
## [441,] "1131873.97181111"   "RAN"      
## [442,] "2225694.00032936"   "RAN"      
## [443,] "2487387.00847225"   "RAN"      
## [444,] "18754.0635224065"   "LOWCARBON"
## [445,] "153238.051869403"   "RAN"      
## [446,] "1900215.30615228"   "RAN"      
## [447,] "1413233.49798306"   "RAN"      
## [448,] "115.430806874736"   "LOWCARBON"
## [449,] "166814.283585256"   "RAN"      
## [450,] "1707727.0507292"    "RAN"      
## [451,] "904634.041106401"   "RAN"      
## [452,] "126458.029982815"   "RAN"      
## [453,] "1192391.15413939"   "RAN"      
## [454,] "1167369.26995129"   "RAN"      
## [455,] "512992.72511025"    "RAN"      
## [456,] NA                   "CRASHED"  
## [457,] "66995.4914433489"   "LOWCARBON"
## [458,] "1500316.82656525"   "RAN"      
## [459,] "778062.924384724"   "RAN"      
## [460,] "19313.4679437398"   "LOWCARBON"
## [461,] "1358187.39438262"   "RAN"      
## [462,] "2141129.84897055"   "RAN"      
## [463,] "7961.46211772983"   "LOWCARBON"
## [464,] "1761228.27665708"   "RAN"      
## [465,] "667103.947334787"   "RAN"      
## [466,] "1516391.6244355"    "RAN"      
## [467,] "2069132.40213798"   "RAN"      
## [468,] "2188.4705746542"    "LOWCARBON"
## [469,] "1156158.26689581"   "RAN"      
## [470,] "3456195.90521345"   "RAN"      
## [471,] "2280381.4764162"    "RAN"      
## [472,] "679642.282835393"   "RAN"      
## [473,] "1.07229600832304"   "LOWCARBON"
## [474,] "2476588.47000487"   "RAN"      
## [475,] "838463.746365504"   "RAN"      
## [476,] NA                   "CRASHED"  
## [477,] "1458452.40164747"   "RAN"      
## [478,] "44044.2655378839"   "LOWCARBON"
## [479,] "651655.234920384"   "RAN"      
## [480,] "1970699.20064016"   "RAN"      
## [481,] "1405.80514674083"   "LOWCARBON"
## [482,] "78.8471668056193"   "LOWCARBON"
## [483,] "14637.1482404126"   "LOWCARBON"
## [484,] "732601.372513896"   "RAN"      
## [485,] "446047.596040638"   "RAN"      
## [486,] "1838908.81647584"   "RAN"      
## [487,] "236713.773033354"   "RAN"      
## [488,] "2479198.89697359"   "RAN"      
## [489,] "1279945.11633116"   "RAN"      
## [490,] "56645.3158685308"   "LOWCARBON"
## [491,] "3255383.23008106"   "RAN"      
## [492,] "597093.74530109"    "RAN"      
## [493,] "1529677.02897881"   "RAN"      
## [494,] "787714.938366425"   "RAN"      
## [495,] "2798.98766749579"   "LOWCARBON"
## [496,] "319.650226473177"   "LOWCARBON"
## [497,] "832705.134925607"   "RAN"      
## [498,] "2919543.18256786"   "RAN"      
## [499,] "513561.801092756"   "RAN"